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ABSTRACT 

Stability and attractor property of free-floating axisymmetric magnetic bubbles in high- 
conductivity plasmas is (tentatively, numerically) demonstrated. The existence of compact non- 
singular axisymmetric magnetic equilibria is proved. Being attractors, the solitary magnetic 
bubbles should exist in nature. 



1. Introduction 

We show that compact non-singular axisym- 
metric equilibrium magnetic field configurations 
are attractors (more precisely intermediate asymp- 
totics) in a weakly dissipative plasma. This means 
that some initial magnetic field configurations 
evolve into solitary magnetic bubbles. 

Because of the attractor property, such free- 
floating solitary magnetic bubbles should exist 
in nature (see Braithwaite 2010 and references 
therein). 

2. Generic and intermediate-asymptotic 
magnetic equilibria 

First consider magnetic equilibria of ideally 
conducting incompressible plasmas. The magnetic 
field B will be in equilibrium if the magnetic j x B 
force is balanced by the pressure gradient: 



B X V X B 



(1) 



It turns out, that due to topological obstacles 
(Arnold 1986), generic magnetic equilibria, that 
is generic solutions of eq. ffl , contain singular cur- 
rent layers (Gruzinov 200yJ). 

If the plasma electric conductivity is finite, but 
still very high, the magnetic diffusivity will dis- 
perse the singular current layers by the reconnec- 
tion process, without dispersing the entire mag- 
netic field configuration. This is the intermediate- 







Fig. 1. — Upper Left: Initial isolines of ij]. Dif- 
ferent colors indicate different signs of ■(/'• Upper 
Right: Singular equilibrium stage. Lower Left: Re- 
connection stage. Lower Right: Non-singular equi- 
librium stage 



^ The main prediction of this paper - generic equihbrium has 
a dense set of singular current layers - has been confirmed 
(tentatively, numerically) 
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asymptotic stage we are interested in. During the 
intermediate-asymptotic stage, the magnetic equi- 
Ubrium configuration should be non-singular. 

Ideally, we want to find all compact non- 
singular magnetic equilibria. We think we have 
a full solution of this problem in two dimensions 
(§3 ). In 2D, the only compact non-singular equi- 
librium is a circular-symmetric magnetic bubble. 

In 3D, there is a class of non-singular equilib- 
rium configurations - compact axisymmetric mag- 
netic bubbles (§4 ). The non-singular axisymmet- 
ric equilibrium bubbles are attractors (§5 ). But 
unlike the 2D case, we do not know if these bub- 
bles are the only non-singular compact attractors. 

We study magnetic equilibria using the relax- 
ation method (e.g. Gruzinov 2009). We put 
an arbitrary compact magnetic field into high- 
viscosity small-magnetic-diffusivity plasma, and 
let the plasma move under the resulting magnetic 
stress: 

B = Vx(vxB)-f?7AB, 
V = -A-^v X V X (B X V X B), ^' 

The first equation (with 77 = 0) says that the 
magnetic field B is frozen into the fluid moving 
with velocity v. Non-zero rj describes magnetic 
diffusivity, A = V^. The second equation says 
that the fluid incompressibly yields to the mag- 
netic force. The equations are solved numerically, 
using the fast Fourier transformation. 

Needless to say, the creeping flow evolution is 
not what happens in nature. But the flnal non- 
singular equilibria that we find should be real. The 
idea is that the route to the equilibrium is irrele- 
vant, we just find generic non-singular equilibria. 

3. Generic and intermediate-asymptotic 
magnetic equilibria in 2D 

Consider 2D magnetic fields, that is magnetic 
fields of the form 



B = {-dyi},d^^p,0), ^^^{t,x,y). (3) 

Take a generic initial magnetic stream function 
■0 vanishing at spatial infinity and let it evolve 
according to eq. ^ with zero magnetic diffusivity. 



?7 = 0. Then each saddle of -0 develops into a 
singular current layer (Gruzinov 2009). 

If we want to have a non-singular equilibrium, 
we must start with the initial ■0 without saddles. 
Then the isolines of "0 are nested closed curves. 
The ideal relaxation of magnetic energy simply 
turns all the nested curves into concentric circles, 
giving the final state -0 = 0(r) where r is the dis- 
tance from some center. 

Now consider slightly non-ideal relaxation, with 
magnetic diffusivity 77 > but still very small. 
Numerical experiments (some 300 different runs) 
show that the generic initial state with compact 
support evolves as follows: 

• Sharp current layers form at saddle points 

• Reconnection occurs at current layers 

• The field splits into a number of circular 
magnetic bubbles. 

• The bubbles have monotonic '0('"): with a 
maximum or minimum at r = and vanish- 
ing -0 at large r. 

• Bubbles of the same polarity (same sign of 
ip) merge, bubbles of the opposite polarity 
repel each other. 

• Ultimately, two bubbles form, one positive 
and one negative, and the distance between 
the two final bubbles is gradually increasing. 

This is illustrated in fig.l. At high numerical 
resolution, when the small magnetic diffusivity can 
be achieved, the different stages of the relaxation 
process are easily distinguished. At first, ideal re- 
laxation quickly moves the field into the minimum- 
energy-for-a-given-topology state. Then reconnec- 
tion occurs, on the time scale which is much longer 
than the ideal relaxation time. After that, an al- 
most steady non-singular equilibrium state is es- 
tablished. 

The 2D summary is: (i) generic equilibrium is 
singular, (ii) the only non-singular compact equi- 
librium configurations are circular bubbles, (iii) 
for arbitrary compact initial field, circular solitary 
magnetic bubbles are the only attractors (more 
precisely, intermediate asymptotics of the weakly 
dissipative MHD). 
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This behavior - formation of two infinitely sep- 
arated circular bubbles ~ seems too simple to be 
true. But in fact its origin is quite clear. In 2D, 
there are infinitely many robust invariants, which 
are constant even when the magnetic energy can 
decrease due to magnetic diffusivity. From eq.©, 

^ J (fr F{i;) ^-vj d^r F"{i:){V^f. (4) 

This is very small, because the rate of dissipative 
magnetic energy damping is given by the integral 
of r/(AV')2. 

The final state is simply the state with min- 
imal energy for the given area per each interval 
dij}. This minimum is given by the state with two 
monotonic circular bubbles, one positive and one 
negative. 

4. Axisymmetric non-singular magnetic 
equilibria 

Now we want to find the weakly-dissipative 3D 
attractors. Unlike the 2D case, we cannot claim 
a full understanding of this problem. This is be- 
cause the runs take longer, and the magnetic fields 
are harder to visualize in 3D. Also, in 2D, we know 
that all saddles become singular current layers un- 
der the ideal relaxation. In 3D, generic ideal equi- 
libria are also singular, with infinitely many singu- 
lar current layers (Gruzinov 2009), but it remains 
unclear where the current layers form. 

However one thing seems to emerge from our 
numerical experiments with the weakly dissipa- 
tive magnetic relaxation in 3D (about 100 differ- 
ent runs). There exist non-singular axisymmet- 
ric compact magnetic configurations, which are 
attractors. In this section we describe the non- 
singular axisymmetric equilibrium bubbles, and 




Fig. 2. — Left: Isolines of from ea.p^ with 
A = 2.74. Right: A = 0.523 



in §5 we demonstrate their attractor and stabil- 
ity properties. 

In cylindrical coordinates {r,9,z), the axisym- 
metric magnetic field can be written as follows 

B = i(-a,^,A,9.V), (5) 
r 

where ip = il}{r^ z), A — A{r, z). Then the current 
is 

V x^ = -{-d,A,-A^il},drA), (6) 
r 

where A* = 9^ — + 5^ . The equilibrium equa- 
tion ([T]) then reads 

l^<,i}drip + AdrA + r'^drP^O, (7) 

A,V9^V + ^9^A + r2a^p = 0, (8) 
dri^d.A - d.il^drA = 0. (9) 

From eq.®, A = A{iIj). Then from eqs.d?] jS]), 

p = piijj), and we get the Grad-Shafranov equation 

A,V + AA' + r^p = 0. (10) 

where ' = ^. 

This equation describes axisymmetric magnetic 
configurations which minimize magnetic energy 

E = j dr dz\[Wij f +A^) (11) 

for a fixed volume between magnetic surfaces 

/[/] = j dr dz rfi^P) = const, V/, (12) 

and a fixed toroidal flux between magnetic sur- 
faces 

J[g] = J dr dz -g{ip)A = const, Wg. (13) 

This variational formulation of the Grad-Shafranov 
equation suggests that non-singular solutions 
should exist for the trivial topology of magnetic 
surfaces. The trivial topology here means the 
absence of separatrix surfaces. The magnetic sur- 
faces should be nested tori with a single maximum 
line. 

The invariant functionals / and J, being lin- 
ear, can be replaced, without loss of generality, by 
some invariant functions of a single variable, say 
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I{a) = - a)], J{a) = J[S{ip - a)] - the vol- 

ume and the toroidal flux distributions functions. 
Thus, we see that the axisymmetric magnetic equi- 
libria are parametrized by two functions of a single 
variable. 

In the 2D case, we only have one invariant func- 
tion /. The final state is also a function of a sin- 
gle variable, ijjir). Thus, it is not surprising that 
the 2D non-singular equilibria are all circular sym- 
metric. In the 3D axisymmetric case, we have two 
constraints, and therefore the shape of the non- 
singular solution cannot be universal. If one keeps 
/ fixed, and changes J, the tori of constant ip will 
incompressibly deform. 

The Grad-Shafranov equation (ITU)) can be eas- 
ily solved numerically (and even analyzed theo- 
retically for some A and p). Perhaps the sim- 
plest non-trivial case (coming from p — —^ijP', 

= iA^2 ^ 1^3) 

- + Vilj^\ip (14) 

Here V — —i} \s the potential of the Shrodinger 
equation, and A is the energy level. The potential 
is confining in the r direction. Confinement along 
z is self-generated by positive ij). It is clear that 
non-singular solutions exist, at least for A slightly 
smaller than the ground state energy of the radial 
part of the hamiltonian ( the operator -d^ + -(?r + 

To confirm the non-singularity, and for the sta- 
bility studies of §5, we have solved eq.([T4]) numer- 
ically, using the relaxation method 

V' ^ C- (V'-l-di(A*V + V'^ -^V)), (15) 

where the time step dt should be small enough, 
and the normalization factor C at each time step 
is chosen so as to keep the maximum value of if) at 
a pre-fixed value. The value of A is given by the 
asymptotic value of C, C = 1 -I- Xdt. The results 
are shown in fig. (2). 

In summary, a family of non-singular axisym- 
metric equilibria parametrized by two functions of 
a single variable exists. 

5. Stability and the attractor property. 

Are the axisymmetric non-singular bubbles sta- 
ble in 3D? Are these bubbles attractors of the 
weakly dissipative MHD? 



5.1. Stability 

Stability of plasma in magnetic fields is the 
most elaborate part of plasma physics (Kadomt- 
sev 1966). The development of the MHD stability 
theory was motivated by the practical problem of 
confining the nuclear burning plasma by external 
magnetic fields. 

In astrophysics, the MHD stability theory has 
found another application. Instead of asking 
which external magnetic configurations will stably 
confine the plasma, one asks what magnetic con- 
figurations can be stably supported by the plasma. 
What are the possible relic magnetic fields in non- 
convective stars, or what kinds of magnetic relics 
can exist in the ISM? 

This problem (formulated by Chandrasekhar 
and Fermi 1953) is still unsolved. One can show 
that certain magnetic field configurations are un- 
stable (Tayler 1982 and references therein). But 
the proof of stability is more difficult. One 
needs to show that the energy functional is non- 
decreasing under an arbitrary virtual displace- 
ment, and this turned out to be hard (although 
Gruzinov 2008 proves stability for some artificial 
magnetic configuration). 

A possible way to check the stability of a given 
magnetic field numerically is by the ideal (meaning 
77 = 0) relaxation method ([2|). If the equilibrium is 
the energy saddle or maximum, rather than min- 
imum, the relaxation will destroy it. We checked 
that the relaxation scheme does catch the classi- 
cal kink and screw instabilities (for definitions, see 
Kadomtsev 1966) 

When we take the equilibrium bubble with 
A = 2.74 (shown in fig. 2), add a random weak 
magnetic field (to seed the possible instabilities), 
and subject it to the 3D relaxation we see just 
the relaxation of the random component (mag- 
netic energy decreases a bit). We also see that 
the noise slightly tilts the bubble (fig. 3). 

Of course we cannot conclude that the bubble 
is stable, especially given the danger of resonant 
surfaces where the field lines are closed. But the 
instabilities, if any, should be mild and globally 
non-destructive. 
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5.2. Attractor property 



6. Summary 



Given that the non-singular axisymmetric bub- 
bles are stable, it is tempting to suggest that they 
behave just like the circular 2D bubbles. This 
is true, except for one major difference. Recall 
that the opposite polarity 2D bubbles repel, and 
the same polarity bubbles merge. In 3D, the 
role of polarity is played by the magnetic helic- 
ity H = J d^rA • B, B = V x A. But unlike the 
2D case, the 3D bubbles always seem to ultimately 
merge, regardless of the signs of helicity and axis 
orientations. 

The difference between the 2D and the 3D cases 
must be due to the fact that in 2D we have in- 
finitely many robust invariants /. In 3D, / and 
J are invariant (and even defined) only under ax- 
isymmetry. The only generic robust invariant in 
3D is, apparently, the helicity. Then, if the total 
helicity is zero, the magnetic field has no reason 
to stay, and indeed it does seem to decay to zero. 
More work is needed here. We can only show what 
we have seen: 



We have shown (tentatively, numerically) that 
non-singular axisymmetric equilibria are weakly 
dissipative attractors. The unsolved problems are: 

1. Do other weakly dissipative (that is non- 
singular) attractors exist? 

2. Do the weakly dissipative attractors form 
(and thus exist) in nature? 
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• Fig. 4. Just like the 2D non-circular bubbles 
relax into circular bubbles, in 3D, the non- 
axisymmetric bubbles relax into axisymmet- 
ric. 

• Fig. 5. Just like the 2D bubbles of the same 
polarity merge into a single circular bubble, 
in 3D, bubbles of the same helicity merge 
into a single axisymmetric bubble, conserv- 
ing helicity but decreasing energy. Fig. 5. 
shows the merger of the bubbles with the 
same helicity and with parallel axefl 

• Fig. 6. The merger of the bubbles with the 
same helicity and with anti-parallel axes. 

• Fig. 7. The 2D bubbles of the opposite polar- 
ity repel. In 3D, the opposite helicity bub- 
bles still merge. If the total helicity is zero, 
apparently a total destruction of the mag- 
netic field occurs. Fig. 7 shows annihilation 
of the opposite helicity bubbles with anti- 
parallel axes. 



^ We define the axis of the bubble to point in the direction 
of the magnetic field on the axis 

This 2-coIumn preprint was prepared with the AAS lATJiiX 
macros v5.2. 
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Fig. 3. — Upper Left: Initial isolines of By in the 
plane y = 0.5. Different colors indicate different 
signs of By Upper Right: Initial B^ in the plane 
z = 0.5. Lower Left: Final By Lower Right: Final 




Fig. 5. — Same as fig. 3. The energy decreases by 
34% . 
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Fig. 7. — Same as fig. 3. The energy decreases by 
factor 13. 
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